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Abstract 



The branching fraction for the decay K ir^eu is measured using in-flight decays de- 
tected with ISTRA+ setup working at the 25 GeV negative secondary beam of the U-70 PS: 
BrK,^ = (5.124 ± 0m9{stat) ± 0.029(norm) ± OmO{syst))%. 

From this value the \ Vus \ element of the CKM matrix is extracted, using previously measured 
f+(t) form factor: \Vus\ = 0.227±0.002. The results are in agreement with recent measurements 
by BNL E865, FNAL KTeV, KLOE . 
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1 Introduction 



The decay K — > euT:^{Kf,z) is known to be one of the best sources of information about Kis 
element of the Cabibbo-Kobayashi-Maskawa (CKM) matrix. The interest in high statistics, 
low systematics measurement of the i^es branching has raised after the paper by BNL E865 
collaboration [T], where 2.5 a increase of the ii'ij branching as compared with PDG02 [2] was 
reported. This result improved the agreement of the measured Vud,Vus,Vub with the unitarity 
condition: 

|Kd|' + |Ks|' + |K6|' = l, 

which was violated by 2.3 a with the old value of Vus- Since then, a set of new measurements of 
the Kf.3 branchings for [3l[4l[5], i^sfSj has appeared, confirming the increase of Vus value. In 
our analysis, we present a new measurement of branching based on statistics of about 2M 
events using new approach, which allows to significantly reduce the systematics uncertainties. 

2 Experimental setup 

The experiment has been performed at the IHEP 70 GeV proton synchrotron U-70. The 
experimental setup ISTRA+ (Fig. 1) has been described in some details elsewhere [7]. 




Figure 1: Elevation view of the ISTRA+ detector. 



The setup is located in the negative unseparated secondary beam. The beam momentum in 
the measurements is ~ 25 GeV with Ap/p ~ 1.5%. The admixture of K~ in the beam is ~ 3%. 
The beam intensity is ~ 3 ■ 10^ per 1.9 sec. U-70 spill. The beam particle defiected by Mi 
is measured by BPCi ^ BPC4 PC's with 1mm wire step, the kaon identification is done by 
Cq -i- C2 threshold C-counters. The 9 meter long vacuumed decay volume is surrounded by 8 
lead glass rings LGi -7- LGg used to veto low energy photons. SP2 is a lead glass calorimeter 
to detect/veto large angle photons. The charged decay products defiected in M2 with 1 Tm 
field integral are measured with the help of PCi PCs — 2 mm step proportional chambers; 
DCi -r- DCs ^ 1 cm cell drift chambers and finally with 2 cm diameter drift tubes DTi 4- DT4. 
Wide aperture threshold Cerenkov counters C3, C^ are filled with He and are not used in the 



2 



present measurements. SPi is a 576-cell lead glass calorimeter, followed by HC — a scintillator- 
iron sampling hadron calorimeter, subdivided into 7 longitudinal sections 7x7 cells each. MH 
is a 11x11 cell scintillating hodoscope, used to improve the time resolution of the tracking 
system, MuH is a 7x7 cell muon hodoscope. 

The trigger is provided by Si S5 scintillation counters, Cq Cerenkoy counters, analog 
sum of amplitudes from the last dinodes of the SPi : T = Si ■ ^2 ■ 6*3 ■ ^4 • Co • • (72 ■ •S's ■ S(S'Pi), 
here 5*4 is a scintillator counter with a hole to suppress beam halo; S'5 is a counter downstream 
the setup at the beam focus; S(S'Pi) — a requirement for the analog sum of amplitudes from 
SPi to be larger than ~700 MeV — a MIP signal. The last requirement serves to suppress the 
K —y iiu decay. 

3 General description of the experimental method 

Our experimental approach to the K^s branching ratio measurement is based on the following 
points: 

1. is the dominant source of electrons in single track decays of charged kaon. Indeed, 
Btk^, ~ 5%;BrK,, ~ 1.5 x IQ-^; Btk^,^ ~ 1.5 x 10"^ fir^^^^^^ ~ 2 x 10"^ The 
contribution from the decay chain K ^u; /i evv corresponds to effective Br^f.^ < 
10~^, because of the long lifetime of muon. The background sources do not exceed fraction 
of % from i^'es and can be easily taken into account. 

2. The number of electrons is obtained from the fit of the Ej-p distribution, where E is the 
energy of the shower, associated with the charged track with momentum p. 

3. The decay K~ — > 7r~7r''(_ft'^2), which is used for the normalization is identified by the 
peak in the distribution over momentum of the charged secondary track in the kaon 
c.m.s (p^™''— distribution) . The number of Kt,2 events is obtained from the fit of this 
distribution. 

4. This method is based on the reconstruction of the beam and decay track only, i.e does 
not require a reconstruction of 7r° both in i^es and K^.^ decays. It uses few selection cuts, 
thus one can hope for a small systematics in this analysis. 

4 Data set and event selection 

During physics run in Winter 2001, 332M events were logged on tapes. This information is 
complemented by about 260M MC events generated with Geant3 [8] Monte Carlo program. 
MC generation includes a realistic description of the setup including decay volume entrance 
windows, track chambers windows, gas, sense wires and cathode structure, Cerenkov counters 
mirrors and gas, shower generation in EM calorimeters, etc. 

The data processing starts with the beam particle reconstruction in BPC\ ^ BPC4, then 
the secondary tracks are looked for in PCi PC3 ; DCi DC^; DTi ^ DT^^ and events with 
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one good negative track are selected. The decay vertex is searched for, and a cut P > 1% on 
the probability of the vertex fit is introduced. 

The decay vertex is selected to locate in the decay volume region 500 < z < 1500 cm, z 
being the coordinate along the beam line. To suppress undecayed particles(beam electrons, in 
particular) a cut on the space angle between beam and secondary track is introduced: A6' > 2 
mrad. 

The next step is to require the total energy in the SPi calorimeter to be above one GeV: 
Espi > 1 GeV. This cut repeats "digitally" the trigger requirement, which is introduced to 
suppress K ^ fiu decays. The matching of the charged track and a shower in SPi is done on 
the basis of the distance r between the track extrapolation to the calorimeter front surface and 
the shower coordinates (r < 5 cm). This cut is used for the electron identification only. 



5 Verification of the method on Monte-Carlo events. 

The cuts described above were applied to the MC-sample which contains a natural mixture of 
reconstructed six largest kaon decays 7r~7r°, n~'7i~^n~, n^e'v, n^fi'i', 7r~7r°7r°), i.e the 

sample includes both signal and main backgrounds. 
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Figure 2: The cumulative distributions over the Figure 3: The cumulative distributions over 



ratio of the energy of the associated calorimeter 
cluster to the momentum of the charged track 
{E/p plot) for four largest background decays 
and Ke3 signal events (MC-events). 



pcms _ momentum of the secondary particle 
in the kaon c.m.s. system, assuming that the 
particle is vr-meson for four largest background 
decays and K^s signal events (MC-events). 



The E/p distribution for these events is presented in Fig. 2. MC shows that the main back- 
ground to E/p is from n~7i^ and vt^tt^tt^. Background is smooth enough to be described 
by A X e~^^'^, signal is described by sum of two Gaussians. Direct test of the fit 



1.006 X Ar*™^ 



gives 

where N^^^^ — is the number of events in the peak of E/p distribu- 



tion of Fig. 2 and N^^^ is the "true" number of Kg^ events in Fig. 2 known from MC. 

The p'^^ distribution for the same events is presented in Fig. 3. For the p^^ the main 
background is cutt^ and /iuti^. Background is smooth enough to be described by 4-th order 
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polynomial. Signal is described by the sum of two Gaussians. 

6 Data Analysis 




Figure 4: The E/p distribution for the real 
data. 




Figure 5: The distributions over the for 
the real data. 



The application of the procedure described above to the real data results in the E/p and 
pcms distributions of Fig. 4, Fig. 5. The fit of the distributions gives: 

Nk,^ = (2.1739 ± 0.0024) x 10^ Nk^^ = (10.2940 ± 0.0053) x 10^ for the data 
Nk,s = (1.2319 ± 0.0013) X 10^ Nk^^ = (6.2758 ± 0.0030) x 10^ for the MC. 

In GeantS version which we are using, the following branchings are assumed: Br^^^ = 
4.82%; Btk^^ = 21.17%. From this we can get: 

BrKes/BrK^^ = 0.2449 ± 0.0004(stat). 

Using PDG06 [9] value BrKn2 = (20.92 ± 0.12)%: 

BrKes = (5.124 ± 0.009(stat) ± 0.029(norm))%. 

In fact, the Kt^2 branching of [9] is obtained by the fit, which has many inputs, including 
BfKesI BrK^2 ratio. That is, it would be more correct to repeat the fit with our new result on 
branching ratio. This is done in [Hj together with averaging over all recent experimental data. 
In present paper we however decided to limit ourself to our own results. 



7 Study of systematics 

The specific feature of our measurements is that the statistical error is much smaller than the 
systematic one. This allows us to study systematics by subdividing our statistics in parts over 
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different variables. Fig. 6 shows the dependence of the measured Nk^.^/Nk^^ ratio versus run 
number and Fig. 7 versus z — the vertex coordinate along the beam line. 
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Figure 6: The measured ratio Nk^-;^/Nk^^ Figure 7: The measured ratio Nk^.;^/Nk^2 

versus run number. versus z coordinate of the vertex. 

The spread of the measured values around average is significantly larger than that ex- 
pected from gaussian statistics. In extracting the systematics error we used an approach 
proposed by PDG |9j when calculating average from different experiments: a scale factor is 
defined s = [x^/iN — 1)]^''^. If we scale up all the errors by this factor, the becomes N-1, as 
required by ideal Gaussian statistics and the error of the average scales up by the same factor. 
The systematic error is then defined as agyst = o'statVs'^ — 1. For example, the scale factor for 
Fig. 6 equals 3.25 and for Fig. 7 it is 3.43. Note, that Fig. 7 demonstrates clear systematics 
in the region of z before the vacuumed decay volume {z < 700 cm) and after it {z > 1300 
cm). We could cut out this regions, reducing systematics related to z, but it does not reduce 
significantly the scale factors for other distributions. That is why we decided not to introduce 
this "a posteriori" cut. 

In this way several more distributions were studied, in particular, over azimuthal and polar 
angles of the secondary track etc. The average scale factor observed is s ~ 3.5. This gives 
estimation for systematic error in BrK^^/BrK^2 of 0.0014 or the systematic error in Btk^^ of 
0.030%. 

A systematics related to a possible admixture of electrons, muons and pions in the beam 
was separately studied. Fig. 8 shows the distribution of events from the E /p peak of Fig. 4 
over momentum (p). The histogram corresponds to the selection of the E/p peak region by a 
simple cut, and the points with errors are the results of the fit of E/p distribution for every 
bin in momentum. The absence of a bump in the region of the beam momentum (~ 25 GeV) 
indicates the absence of this type of background. Indeed, beam electrons are suppressed by two 
C-counters in the beam and by the > 2 mrad cut. 

It is easy to show that the amount of electrons from the decay chain tt — > /i — > e of the 
beam pion is negligible as compared with the decay chain K fi ^ e, which is correctly 
reproduced in our MC: the number of pions which may pass C*-counters "Veto" is at most 
0.5%, i.e it is 1/6 of kaons, the lifetime factor (7cr) is 7.5, then taking into account the K ^ fiu 
branching of .63 we get factor of 30 in favour of kaon decays. 
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Figure 8: The measured momentum distribution of the tracks, identified as electrons. 
Histogram corresponds to the "rough" identification by the selection 0.6 < E/p, the 
points with errors — results of the fits of the E/p distributions for each bin in p. 

The effects on Nk^^/Nk^2 from the cuts variation (vertex fit probability, > 2 mrad, Espi > 
1 GeV) and different parametrization of the signal and background of Fig. 2 - Fig. 5 are less 
then estimation of systematic error from previous section. 
Summing up all the systematics observed leads to our final result: 

BrK^jBvK^^ = 0.2449 ± 0.0004(stat) ± O.OOU{syst) , 

Btk,^ = (5.124 ± 0.009(stat) ± 0.029(norm) ± 0.030{syst))% . 
The comparison of our results with the E865 pl shows reasonable agreement. 



8 Extraction of \Vus\ 

In the Standard Model the Born-level matrix element for the tt^I^u decay modes is: 

here ^[/r^" " {PK+Pn)a + /5^"" ■ {pk - P.)a] =< Asiail - lMK+ > ; - form 

factors, which depend on t = {pk — P-kY = {Pi + PvY — the square of the four momentum 
transfer to the leptons. 

The term in the vector part, proportional to /_ is reduced (using the Dirac equation) to an 
effective scalar term, proportional to mi. That is why in case of Kes decay one can neglect the 
term proportional to /_ . 

The Kf^, decay rate can be expressed as: 

r(i^e1) = = ^,Ml,\Vus?\Um'lK^SEw{l + 5sU2 + 5X? . 
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Here Sew = 1.0232 ± 0.0003 is the short-distance radiative correction [10] ; 
Ssu2 = (2.31 ±0.22)% takes into account the difference between ff'^'iO) = /+(0) and (0) 
[TTl [T4] : 5^ = (0.03 ± 0.1)% is the long distance radiative correction for K^^, for the fully in- 
clusive A'e3(7) decay The is the dimensionless decay phase space integral |12j: 

Ik^ = I dtj^X'/\U{t)/U{0)r . 

Where A = (M| - t - Mlf - AtM^ 

To extract \Vus\, we take Brj^±^ from the present experiment, t{K^) = (12.385 ± 
0.025) nsec — the charged kaon life-time from PDG06 [9] and calculate from our mea- 
surement of /+(t), where the quadratic non-linearity was observed for the first time |13| : 

Fk+ = 0.15912 ± 0.00084(stat) ± 0.00lU{syst) . 

The systematic error reflects the difference between the quadratic and linear flt of the f+(t). 
Putting everything together we get: 

|V;j+(0)| = 0.2186 ± O.OOOQb^ ± 0.0012t;, . 

And finally: 

\Vus\ = 0.2275 ± 0.0009b^ ± 0.0022^;, . 
If theoretical value /+(0) = 0.961 ± 0.008 pjj is used. 

9 Summary and conclusions 

The decay has been studied using in-flight decays of 25 GeV K~ , detected by ISTRA+ 
magnetic spectrometer. Due to the high statistics, adequate resolution of the detector and good 
sensitivity over all the Dalitz plot space, the errors are signiflcantly reduced as compared with 
the previous measurements. The Kf,^ branching is measured to be: 

BrK^jBr-K^^ = 0.2449 ± 0.0004(stot) ± O.OQU{syst). 

Btk,, = (5.124 ± 0.009(stat) ± 0.029(norm) ± 0. 030 {syst))%. 

From that we obtain: 

\Vusf+{0)\ = 0.2186 ± 0.0009b^ ± 0.0012^/,. 

Which leads to: 

iKsl = 0.2275 ± 0.0009Br ± 0.0022^;, 

if the theoretical value for /+(0) is substituted. Our result on \Vus\ is in reasonable agreement 
with that from charged [ij and neutral [3], [5], [6] kaon decays. 

The work is partially supported by the RFBR grants N07-02-00957, N07-02-16065 and by 
the grant of the Russian Science Support Foundation. 
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